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EXECUTIVE  SUMMARY 


OBJECTIVE 

Assess  the  accuracy  of  the  Sudden  Ionospheric  Disturbance  Model  (SIDGRID) 

and  present  current  computer  implementations  of  the  algorithm. 

RESULTS 

1.  The  NOSC  solar  X-ray/HF  lowest  observable  frequency  (  LOF )  algorithms  provide 
an  accurate  assessment  of  the  change  in  HF  LOF  for  solar  X-ray  bursts  greater  than 
1  x  10— ^  erg/cm^  x  s  and  are  valid  over  the  entire  sunlit  hemisphere. 

2.  The  SIDGRID  model  has  the  ability  to  display  the  disturbed  lowest  usable  frequency 
(LUF)  for  any  path  on  the  geographical  sunlit  hemisphere  at  any  date  and  time. 

3.  For  distances  less  than  1500  km,  the  short-range  prediction  error  was  ±0.5  MH/. 

4.  For  distances  greater  than  1500  km.  the  long-range  prediction  error  was  less  than  or 
equal  to  ±2.0  MHz  in  95%  of  the  tests. 

RECOMMENDATIONS 

1.  Modify  the  SIDGRID  model  to  include  a  means  to  adjust  the  resultant  LUFs  for 
system  parameters  such  as  transmitter  power,  antenna  gains,  required  signal-to- noise 
ratio,  and  background  radio  noise. 

2.  Modify  the  SIDGRID  model  to  improve  the  transition  from  disturbed  to  quiet  time 
LUF  prediction. 

3.  Develop  the  capability  to  calculate  the  D-region  absorption  arising  from  a  short-wave 
fade  (SWF). 


4.  Improve  the  model  by  eliminating  the  discontinuity  at  1500  km.  and  make  the  model 
continuous  as  a  function  of  range. 
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1.0  INTRODUCTION 


The  effective  operation  of  long-distance  high-frequency  (  HF)  communications 
systems  has  increased  in  proportion  to  the  ability  to  predict  variations  in  the  ionosphere. 
These  variations  are  affected  in  a  complex  manner  by  solar  activity  and  seasonal  and 
diurnal  changes,  as  well  as  by  latitude  and  longitude.  Such  a  predictive  capability  has 
permitted  communicators  to  optimize  frequencies,  antennas,  and  other  circuit  parameters. 

Initially,  manual  methods  were  developed  for  analyzing  ionospheric  variations  on 
HF  circuits  of  short,  intermediate,  and  long  distances.  Because  the  manual  methods  were 
laborious  and  time-consuming,  various  organizations  developed  computer  programs  to 
analyze  HF  circuit  performance.  A  commonly  predicted  parameter  in  these  programs  is 
the  maximum  usable  frequency  (MUF).  The  MUF  is  the  highest  frequency  that  can  be 
propagated  by  ionospheric  refraction  between  points  at  a  given  time.  Another  commonly 
predicted  parameter  is  the  lowest  usable  frequency  ( LUF ).  The  LUF  is  the  lowest  usable 
frequency  propagated,  and  is  determined  by  the  amount  of  D-region  absorption.  The 
LUF  over  any  circuit  path  is  established  as  a  function  of  total  path  absorption  with 
respect  to  such  HF  system  parameters  as  transmitted  power,  signal-to- noise  ratio,  and 
antenna  gains. 

More  recently,  the  Naval  Ocean  Systems  Center  (NOSC)  has  developed  a  series 
of  ionospheric  prediction  programs  that  will  run  on  portable  microcomputers. 
MINIMUF-3.  MINIMUF-3.5.  and  MINIMUF-85  are  examples  of  a  series  of  NOSC- 
developed  MUF  models  (ref  1-3). 

This  report  will  describe  results  of  the  uncertainty  assessment  of  the  SIDGR1D 
model  used  to  make  LUF  predictions  for  a  disturbed  ionosphere.  The  accuracy  of  both 
short-path  and  long-path  predictions  of  the  SIDGRID  model  will  be  discussed.  The 
FORTRAN  implementation  of  the  SIDGRID  model,  routine  DLOF.  is  listed  in 
appendix  C  and  the  BASIC  SIDGRID  algorithm  is  listed  in  appendix  D. 

The  SIDGRID  model  was  compared  to  198  observed  sounder  data  points. 

Results  from  these  tests  indicate  that  the  mean  long-range  prediction  error  was  0.24  MHz 
and  the  error  between  predicted  and  observed  data  was  negligible  in  29%  of  the  tests. 

The  probability  that  the  prediction  error  was  less  than  or  equal  to  ±  1 .0  MHz  was  83% 
and  the  probability  of  error  less  than  or  equal  to  ±2.0  MHz  was  95%.  For  distances  less 
than  1500  km.  the  short-range  prediction  error  was  ±0.5  MHz. 


2.0  COMMUNICATIONS  DISTURBANCES  CAUSED 
BY  SOLAR  FLARE  ACTIVITY 

2.1  DISTURBANCE  PHENOMENA 

Any  system  that  is  ionosphere-dependent  is  also  dependent  on  solar  activity  and  is 
subject  to  transient  variations  in  solar  radiation.  There  arc  two  states  of  activity  on  the 
Sun  that  are  of  interest  —  the  quiet  state  and  the  disturbed  state.  There  are  three  broad 
categories  of  radiation  that  are  of  interest  —  X-ray.  ultraviolet  (  UV).  and  particle 
emission.  It  is  the  abrupt  change  from  the  quiet  to  the  disturbed  state,  which  occurs 
during  a  flare  on  the  Sun's  surface,  that  causes  a  transient  change  in  radiation  levels  and 
consequently  causes  abrupt  changes  in  the  Earth's  ionosphere. 

The  quiet  Sun  radiates  energy  over  nearly  the  entire  electromagnetic  spectrum. 

Each  portion  of  the  spectrum  contributes  to  or  affects  the  ionosphere  in  some  way. 

Certain  energy  bands  within  each  type  of  radiation  provide  the  primary  constituent  for 
ionization  at  different  levels  within  the  Earth's  atmosphere.  This  stratification  creates 
distinct  layers,  which,  together,  form  the  ionosphere.  The  lowest,  or  D-region  (50-90  km 
above  the  Earth),  is  believed  to  be  formed  by  X-rays  in  the  1-10-A  spectrum  and  by  Lyman 
a.  UV,  and  particle  radiation.  The  E-region  (90-120  km)  is  created  by  X-rays  in  the 
10-100- A  range  and  UV  radiation  between  100  and  1500  A.  The  F-region  (above 
150  km)  is  formed  by  solar  UV  irradiation  between  170  and  91 1  A. 

Solar  flares  are  the  source  of  most  of  the  variable  radiation  that  has  important 
effects  on  the  terrestrial  and  interplanetary  environment.  A  solar  flare  is  a  localized 
chromospheric  brightening  on  the  sun  and  is  sometimes  accompanied  by  X-ray  bursts, 
energetic  nuclear  particle  emission,  and  radio  frequency  emissions.  The  size  of  the  solar 
eruption  will  dictate  the  degree  to  which  the  various  ionospheric  regions  are  changed  and 
the  subsequent  disruption  of  any  ionosphere-dependent  service. 

During  periods  of  solar  disruption  (flares,  etc. ).  the  energy  levels  of  these  three 
radiation  forms  abruptly  change,  which  alters  the  structure  of  the  ionosphere,  changing  the 
overall  transfer  function  of  the  service  propagating  medium.  These  transient  energy 
variations  have  w'ide  and  varying  effects  on  ionosphere-dependent  systems.  There  are 
two  categories  of  disturbance  phenomena:  ( 1 )  short-term  effects,  which  occur  in  close 
time  correlation  with  the  solar  eruption  and  last  for  several  hours  after  the  event:  and  (2) 
long-term  effects,  which  commence  24-36  hours  after  the  event  and  last  several  days. 
Short-term  effects  are  usually  the  result  of  enhanced  X-ray.  UV.  and  high-energy  particle 
radiation,  while  long-term  effects  are  caused  by  low-energy  particle  bombardment. 

Figure  1  shows  the  time  event  relationship  of  these  phenomena. 

The  phenomena  shown  in  figure  1  are  summarized  below  ( re f  4  ) . 

SHORT-WAVE  FADE  (SWF).  A  short  wave  fade  is  a  decrease  or  loss  in 
signal  strength  of  3-32-MHz  radio  waves,  caused  by  increased  absorption  due  to  1-10  A 
X-ray  enhancement  in  the  D-region  and  10- 100- A  X-ray  enhancement  in  the  lower 
E-region.  Onset  is  in  close  coincidence  with  the  solar  eruption  and  the  effects  can  last 
from  several  minutes  to  several  hours,  depending  on  the  magnitude  of  the  X  ra\  burst  and 
the  HF  frequency  of  operation. 
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Figure  1  Disturbance  phenomena  i relative  timesi 


SUDDEN  PHASE  ANOMALY  (SPA).  An  SPA  occurs  when  the  lom/atu  n  o( 
the  D-region  is  enhanced  by  1  — 10- A  X-rays  lowering  the  effect  i\ e  rellection  height  I  his 
effect  changes  the  phase  of  a  sky-wave  radio  signal  and  is  predominantly  realized  in  the 
3-30-kHz  VLF  frequency  region.  SPAs  occur  in  close  coincidence  with  the  solar  event, 
and  effects  can  be  noticed  1  to  2  minutes  prior  to.  and  concurrently  with,  the  HE  SWL 

SUDDEN  ENHANCEMENT  OF  ATMOSPHERICS  (SEA).  SKA  is 
another  manifestation  of  D-region  enhancement  by  solar  X-rays  and  is  denoted  by  an 
increase  in  low-frequency  atmospheric  noise  in  the  VLF  region  (that  is.  thunderstorms  at 
around  30  kHz). 

SUDDEN  COSMIC  NOISE  ABSORPTION  (SCNA).  The  SCNA  is  a 
sudden  decrease  in  signal  strength  of  galactic  cosmic  noise  received  at  frequencies  high 
enough  to  penetrate  the  ionosphere,  yet  low  enough  to  encounter  measurable  absorption. 
This  phenomenon  is  characteristic  of  the  18-30-MHz  frequency  spectrum. 

SUDDEN  ENHANCEMENT  OF  SIGNAL  (SES).  An  SES  is  a  sudden 
increase  in  VLF  signal  strength  caused  by  D-region  enhanced  ionization.  Enhancement 
of  the  F- region  by  UV  radiation  is  possible,  although  it  is  generally  masked  by  D-region 
ionization. 


SUDDEN  FREQUENCY  DEVIATION  (SFD).  An  SFD  is  an  impulse 
deviation  in  the  frequency  of  2-20-MHz  radio  waves  reflected  from  the  F-region  of  the 
ionosphere.  This  phenomenon  is  caused  by  ionization  enhancements  in  the  upper  D-.  E-. 
and  Fj-regions  caused  by  radiation  above  100  A  and  impulsive  X-rays  in  the  4-  20- A 
wavelengths.  It  is  characterized  by  a  sharp  peak  deviation  between  0. 1  and  61  Hz  (the 
average  is  approximately  1/2  Hz)  and  a  decay  back  to  zero  deviation.  The  duration  is 
approximately  15  minutes  and  is  coincident  with  event  onset. 

SUDDEN  IONOSPHERIC  DISTURBANCE  (SID).  SID  is  a  general  term 
denoting  the  effects  of  enhanced  ionization  in  the  D-region  and  lower  E-region  by  1-10-A 
X-rays  over  the  entire  sunlit  hemisphere  of  the  Earth.  The  SWF.  SPA.  SFD.  SFS.  and 
SCNA  are  all  different  manifestations  of  the  SID.  depending  on  the  portion  of  the 
ionosphere  involved  or  the  frequency  of  interest.  SIDs  are  short  term  disturbance 
phenomena  occurring  in  close  coincidence  with  the  solar  event. 

IONOSPHERIC  STORMS.  The  ionospheric  storm  is  a  long-term  phenomenon 
resulting  from  a  disruption  of  the  Earth's  magnetic  field  by  particles  (electrons  and 
protons)  ejected  from  a  solar  Hare,  and  traveling  at  a  speed  of  1000-2000  km  min.  The 
onset  can  be  anywhere  from  8  to  48  hours  after  a  solar  event  and  primarily  affects  the 
upper  F-region.  It  is  characterized  by  degraded  HF  radio  propagation,  high- latitude 
radio  blackout,  and  sporadic  E.  The  effects  can  last  anywhere  from  8  to  72  hours.  At 
high  latitudes,  ionospheric  storms  create  such  phenomena  as  auroral  displays  and  polar 
cap  absorption  (PCA)  events. 

POLAR  CAP  ABSORPTION  (PCA).  The  PCA  is  a  result  of  high-energy 
particle  influx  from  a  solar  event  and  is  evidenced  by  enhanced  absorption  of  radio  waves 
in  the  polar  regions.  The  PCA  starts  from  1-4  hours  after  the  solar  event  and  can  last 
from  1-3  days. 

2.2  SWF  FREQUENCY  DEGRADATION 

The  characteristics  of  SWFs  and  their  effects  on  HF  systems  have  been  discussed 
in  references  5-8.  For  an  HF  user,  the  short-wave  fade  can  produce  short-term 
degradation  that  varies  from  minor  signal  loss  to  a  total  frequency  blackout  lasting  up  to 
several  hours.  The  effects  of  the  SWF  are  abrupt  and  in  close  time  correlation  with  the 
actual  solar  event.  Typically,  the  duration  of  the  SWF  is  1-4  hours.  The  mechanisms  of 
a  short-wave  fade  are  reasonably  well  understood.  Sky-wave  propagation  via  the  F- 
region  passes  through  the  lower  D-region  and  is  subject  to  absorption.  During  a  solar 
disturbance.  X-ray  emissions  irradiate  the  ionosphere,  the  depth  of  penetration  depending 
on  the  spectrum  of  the  X-ray  flux,  the  solar  zenith  angle  on  the  transmission  path  and 
the  magnitude  of  the  X-ray  burst.  If  the  penetration  is  of  sufficient  magnitude,  the 
D-region  ionization  is  enhanced,  increasing  absorption  in  this  region  of  all  signals  passing 
through  it.  The  immediate  effect  to  the  communicator  is  signal  degradation  or  loss  at  the 
low  end  of  the  usable  spectrum,  and  the  LOF  markedly  increases  1  to  3  minutes  after  the 
onset  of  the  event. 


The  SWF  event  is  a  daylight  phenomenon,  the  intensity  of  which,  on  a  given 
transmission  path,  depends  on  solar  zenith  angle  and  the  peak  magnitude  and  duration  of 
the  solar  X-ray  burst.  The  peak  frequency  penetration,  or  LOF  increase,  and  the 
duration  of  the  spectrum  outace  can  be  directlv  correlated  to  the  shape  and  magnitude  of 
1-8-A  X  -ray  flux.  The  penetration  of  the  LOF  into  the  usable  spectrum  is  a  function  of 
the  magnitude  of  the  X-ray  burst  and  the  solar  position  along  the  transmission  path.  A 
major  event  (>  10  “  *  erg  cm-  ■  s)  will  typically  penetrate  through  about  3  4  of  the 
usable  spectrum.  For  a  minor  event  ( I0_l  erg/cm-  •  s>  X-ray  flux  >  10“-  erg/cm^  ■  s). 
penetration  will  occur  through  half  the  usable  spectrum. 

In  the  HF  spectrum,  the  D- region  acts  as  a  variable  attenuator,  with  the  total 
amount  of  absorption  depending  on  (  1  >  the  number  of  times  the  signal  passes  through  the 
D-region  over  a  circuit  path  and  (2  )  the  launch  angle  of  the  signal.  During  daylight 
undisturbed  conditions,  the  lower  limit  of  the  propagating  spectrum  |  i. e. .  the  lowest 
observed  frequency  (  LOF)|  is  determined  by  the  amount  of  D-region  absorption.  In  the 
strictest  sense,  the  LUF  over  any  circuit  path  is  also  a  function  of  such  HF  system 
parameters  as  transmitted  power,  signal- to- noise  ratio,  and  antenna  gain.  For  this 
reason,  measured  values  of  LOF  will  vary  slightly  for  different  paths  and  HF  systems 
during  undisturbed  ionospheric  conditions.  When  a  solar  Hare  occurs,  the  increased 
absorption  in  the  D-region  becomes  the  dominant  factor  in  determining  the  LOF.  causing 
it  to  rise.  The  amount  it  rises  is  proportional  to  the  peak  value  of  the  X-ray  enhancement 
and  the  solar  zenith  angle.  A  valid  indicator  of  the  impact  of  an  SWF  on  a  system's 
operation  is  the  upward  shift  of  the  LOF. 

For  a  given  HF  circuit,  there  exists  a  band  of  frequencies  that  will  support 
propagation.  The  upper  limit  is  bounded  by  the  maximum  observed  frequency  (  MOF). 
which  depends  on  the  electron  frequency  of  the  E-  or  F-rcgion.  whichev  er  is  supporting 
propagation  at  that  time.  The  lower  boundary  ,  called  the  lowest  observed  frequency 
( LOF).  is  controlled  by  the  amount  of  nondeviative  absorption  in  the  ionospheric  D- 
region.  The  band  of  frequencies  between  LOF  and  MOF  is  called  the  existing  usable 
spectrum.  The  amount  of  absorption  incurred  as  the  signal  passes  through  the  D-region  is 
controlled  by  solar  radiation.  As  the  sun  rises,  the  electron  density  in  the  D-region 
increases,  causing  the  HF  absorption  to  increase  proportionally  to  the  inverse  square  of 
the  propagating  frequency.  Therefore,  during  the  day.  absorption  in  the  lower  HF 
spectrum  increases  and  the  LOF  rises  as  shown  in  figure  2  a.  The  rise  of  the  day  light 
EOF  follows  a  secant  solar  zenith  angle  (\)  function. 

During  solar  disturbances,  the  levels  of  solar  X-rays  rise,  causing  D-region 
absorption  to  increase.  This  has  the  effect  of  causing  the  LOF  to  rise  proportionally  to 
the  increase  in  absorption.  When  the  LOF  meets  the  MOF.  the  circuit  becomes 
completely  blacked  out.  The  phenomenon  of  the  abnormal  increase  in  LOF  caused  by 
the  short-wave  fade  (  SWF )  is  illustrated  in  figure  2  b.  Also  shown  is  the  I  A  X  ray 
burst  measured  by  the  SOLRAD  9  satellite.  The  close  correlation  between  changes  in 
the  HF  LOF  and  1  8-A  X-ray  flux  can  be  easily  seen. 
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The  solar  zenith  angle  is  one  of  the  major  factors  in  determining  the  degree  to 
which  ionization  takes  place  in  the  lower  atmosphere.  The  increase  in  angle  of  incidence 
of  the  solar  radiation  imposes  a  greater  path  length  for  a  given  penetration  of  the  X-ray 
flux  into  the  ionosphere.  A  propagation  path  with  reflection  points  near  the  terminator 
(day-night  transition)  will  exhibit  much  less  signal  absorption  during  a  solar  X-ray  event 
than  will  a  similar  path  with  reflection  points  near  the  subsolar  location.  This  is  because 
the  X-ray  flux  penetrating  near  the  terminator  must  pass  through  considerably  more 
ionosphere  (because  of  the  steep  angle  of  incidence)  before  reaching  the  D- region  where 
the  HF  absorbing  ionization  must  be  generated.  A  large  portion  of  the  X-ray  energy  is 
therefore  lost  before  reaching  the  D- region.  At  the  subsolar  point,  the  shortest  path  to 
the  D-region  is  present  and  the  greatest  quantity  of  X-ray  penetration  is  achieved. 


3.0  SUDDEN  IONOSPHERIC  DISTURBANCE  MODEL  (  SIDGRID) 

SIDGRID  is  an  extension  of  the  development  of  a  simple  LOF  model  that  could 
be  varied  as  a  function  of  real-time  solar  X-ray  flux  measurements  during  solar  flares. 
This  model  was  the  first  truly  real-time,  sudden  ionospheric  disturbance  ( SID)  assess¬ 
ment  capability  for  a  system  like  PROPHET  (ref  0  ).  The  close  correlation  between  the 
concurrent  rise  in  solar  1-A  X-ray  flux  and  HF  LOF  led  to  the  development  of  a  simple 
first-order  algorithm  relating  the  two  as  a  function  of  solar  zenith  angle. 

During  development  of  the  algorithm,  it  became  evident  that  variation  in  the  LOF 
versus  X-ray  flux  dependency  on  transmission  path  length  was  such  that  a  single,  range- 
independent  model  would  not  satisfy  all  requirements.  Analysis  showed  that  improved 
results  could  be  obtained  if  the  transmission  paths  were  separated  into  two  categories  - 
short-range  for  anything  less  than  1500  km  and  long-range  for  everything  greater.  Two 
different  algorithms  were  developed,  one  for  each  range  category.  A  brief  theoretical 
discussion  of  the  formulae  used  to  derive  the  algorithms  is  given  in  appendix  A. 

3.1  METHOD  OF  CALCULATION 

The  long-range  LL'F  X-ray  flux  model  is  an  empirically  derived  relationship.  The  first 
step  was  to  plot  the  l-S- A  X-ray  flux,  as  measured  by  the  SOI. RAD  9  satellite,  versus  the 
I.OF  during  the  onset  of  two  SWFs  observed  over  the  Hawaii-to-California  oblique  sounder 
path  on  30  and  31  October  1968.  These  two  events  were  chosen  because  of  similar  rise  and 
peak  characteristics  of  the  1-8- A  flux,  as  shown  in  figure  3.  The  difference  in  LOF  noted 
between  the  two  curves  arises  from  the  different  solar  zenith  angle  (x)  when  the  events 
occurred.  The  next  step  was  to  derive  a  function  to  account  for  x-  By  “trial  and  error” 
techniques,  the  LOF  at  a  given  time  was  multiplied  bv  a  factor  of  (1.0  +  sec;x  10).  which 
moved  the  curves  together.  By  curve-fitting  techniques,  the  relationship  between  true  l.l  F 
corrected  for  zenith  angle  ( F')  and  the  I -8- A  X-ray  flux  was  determined,  to  a  first  order  of 
approximation,  to  be: 


Flux  (  I  8  -  A  )  (0  01 038  )  |  F'  15). 
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Figure  3.  X- ray  flux  vs  HF  LOF. 


Further  investigation  expanded  this  equation  to  provide  a  better  fit.  the  final 
equation  being: 


Flux  -  0.01038  (f|  15.0)  +  0.003  sin  [0.8491  (f ,  15.6)]  =  0 


where 


f,  =  LOF 


[l  +  sec"  x  j 

L  io  J  • 


To  find  the  LUF,  equations  ( 1 )  and  (2 )  are  solved  by  Newton's  method,  using  the  given 
values  of  x  and  LUF.  Newton's  method  refines  an  initial  guess  xq  of  a  solution  of  a 
general  nonlinear  equation  fix)  =  0.  It  takes  the  following  form: 


<i  + 1  =  *i-j — 

arfu’ 


(i  -  0.  1.  2.  ...  ) 


where  Xg  =  (Flux  0. 1038  +  150)  ( 10  +  sec^x) 

In  equation  ( 1 ).  x  =  LUF.  To  use  Newton's  method,  we  evaluate  d/dx  fi  x)  and  get: 

f.  ',0.01038  -  0.0025473  cos|0.8491  (f,  15.6 )|! 


iAWl;W.>.w 


To  find  the  subsolar  point,  the  Greenwich  hour  angle  (GHA)  and  solar  declination 
are  used.  The  solar  zenith  angle  ( x )  is  computed  from 

X  —  cos-'  | sin  (a)  sin  (8)  +  cos  (a)  cos  (8)  cos  (y  —  fi)\  (5  ) 

where 

a  =  latitude  of  grid  point 
p  =  longitude  of  grid  point 
8  =  latitude  of  subsolar  point  =  declination 
y  =  longitude  of  subsolar  point  =  15  GHA. 

The  short-range  algorithm  was  developed  by  the  same  techniques  discussed  above 
for  the  long-range  algorithm,  but  with  a  different  data  set.  An  extensive  search  did 
produce  a  period  in  1971  when  oblique  sounder  strip  chart  recordings  were  made  of  the 
now-defunct  Davis.  California,  sounding  transmitter.  Enough  information  was  assembled 
to  derive  the  desired  algorithm. 

Davis  to  La  Posta.  California,  is  an  803 -km.  north- to- south  path.  Between  2 
July  and  20  November  1970.  10  HF  SWF  disturbances  were  recorded.  Table  1 
identifies  the  disturbances  analyzed. 


Table  1.  1971  SWF  disturbances  -  Davis  to  La  Posta.  CA  . 


Event 

Date 

Time 

1-8- A  X-ray  Peak* 

1 

2  Jul  1971 

1845 

2.2  X  I0~2  erg/ cm2  •  s 

2 

7  Jul  1971 

1507 

6.3  X  10  2  erg  cur  •  s 

3 

7  Jul  1971 

1603 

8.9  X  10  2  erg  cur  •  s 

4 

7  Jul  1971 

1925 

2.8  X  10"  2  erg'em2  ■  s 

5 

7  Jul  1971 

1947 

5.3  X  10"  2  erg/cm2  •  s 

6 

10  Jul  1971 

1546 

1 .3  X  10" 2  erg/cm2  •  s 

7 

1  2  Sep  1 97  1 

1624 

4.4  X  10"-  erg/ cm-  •  s 

8 

10  Nov  1971 

0000 

4  X  io  2  erg/ cm2  •  s 

9 

12  Nov  1971 

1833 

2  X  10” 2  erg/ cm2  •  s 

10 

20  Nov  197! 

2049 

1.8  X  10" 2  erg/cm2  ■  s 

♦Derived  from  SOLR.AD  9  satellite  1-8- A  X-ray  data. 


The  new  algorithm  was  tested  against  all  available  short-range  SWF  data,  and  it 
became  obvious  that  the  basic  SIDGRID  model  would  require  modification  to  provide  the 
dual- algorithm  capability.  On  completion  of  the  necessary  modifications,  the  over-all 
concept  was  tested  by  use  of  an  SWF  disturbance  that  was  measured  simultaneously  over 
long-range  and  short-range  paths.  Event  8  in  table  1  is  such  an  event. 

From  theoretical  considerations,  the  form  of  the  formula  can  be  given  by 


where 

F  =  LUF  (MHz) 
ip  =  solar  zenith  angle 
Fx  =  X-ray  flux. 

The  contributions  of  energy  other  than  solar  X-rays  are  assumed  to  be  a  constant  effect. 
The  values  of  n,  m.  and  a  are  to  be  found  by  use  of  the  observations.  A  theoretical  value 
for  n  can  be  found  by  noting  the  following  simplified  model  for  absorption.  Electron 
production  rates,  q.  in  the  D-region  of  the  ionosphere  can  be  given  by 

q  =  Cj  Fx  cos'71  V'.  ( 7 ) 

and  q  is  related  to  the  electron  density  at  equilibrium  by 

q  =  aN*.  (8  ) 

The  absorption  at  a  point  in  the  ionosphere  at  which  the  electron  collision  frequency  is  v 
can  be  calculated  from 

A  =  C2N  v/f^.  (9  ) 

It  is  assumed  that  total  absorption  is  proportional  to  some  weighted  average  electron 
density,  N,  which  is  related  to  the  1-A  flux  by  equation  8.  By  combining  equations 
6-9,  it  is  found  that 

A^  =  cFx  cosm  ip  ^  .  (10) 

where  A  is  an  absorption  cutoff  level  determined  by  noise  and  signal-processing  con¬ 
siderations.  The  value  of  n  =  4  has  been  found  to  fit  the  data  quite  well.  To  find  m.  we 
reduce  equation  6  to  a  linear  equation  in  cosm  i p  by  taking  logarithms 

In  ( fxf4)  =  m  •  In  ( cos  ip  )  +  In  (  a).  (11) 

Solving  for  a  and  m,  we  obtain  a  formula  of  the  form  of  equation  10.  in  which  the 
constant  1.0191  x  10  3  is  the  absorption  cutoff  level: 

(1.091  x  10-3)2  =  r4(Fxcos3  1^).  (12) 

It  was  found  that  when  1  -8 A  X-ray  flux  was  below  1.85  X  10"  3  erg/cm3  •  s.  a  "quiet" 
or  undisturbed  approximation  could  be  used: 

f4  =  1785  cos3  iji.  (13) 

This  formula  was  derived  from  pre-event  undisturbed  data  shown  in  table  2. 


3.2  DESCRIPTION  OF  DATABASE 


The  long-range  algorithm  was  derived  by  using  data  from  two  SWF  events  over 
two  HF  transmission  paths.  To  ascertain  how  valid  this  algorithm  is  for  other  paths  and 
different  times  and  seasons,  the  following  four- step  process  was  performed. 

Step  I.  From  the  NOSC  La  Posta  historical  data  files  of  HF  oblique  sounder 
data.  HF  short-wave  fades  were  selected  for  the  period  between  1969  and  1972.  Two 
types  of  SWF  disturbance  data  were  chosen:  ( 1 )  multiple  SWFs  observed  over  a  single 
transmission  path  in  the  same  24-hour  period:  and  (2)  a  single  or  double  SWF  event 
observed  over  more  than  one  path  in  the  same  24 -hour  period.  At  this  stage,  only  the 
facts  that  an  SWF  was  evident  on  the  HF  sounder  recorder  and  that  the  data  were 
suitable  for  further  analysis  were  noted.  No  frequency' time  data  were  recorded  at 
this  time. 

Step  2.  Once  the  date  and  time  of  an  SWF  was  chosen,  a  1-8-A  X-ray  flux 
versus  time  plot  for  the  same  period  was  prepared  from  the  NRL  SOLRAD  9  satellite 
solar  radiation  records. 

Step  3.  After  the  1-8-A  X  -ray  flux  levels  were  known  as  a  function  of  time,  the 
proper  LOF  paths  were  selected  and  the  HF  transmissions  w'ere  plotted.  Under  the 
assumption  that  the  highest  LOF  indicated  along  the  plotted  HF  path  will  control  the 
overall  path  LOF.  a  plot  of  predicted  LUF  versus  time  was  prepared. 

Step  4.  Once  the  predicted  LUF/time  plot  was  completed,  the  HF  sounder 
records  were  analyzed  and  the  observed  LOF  versus  time  was  plotted  for  the  same 
events  In  this  manner,  the  SIDGRID  model  was  tested  in  a  quasipredictive  mode,  with 
only  time  and  X-ray  flux  levels  being  used  to  vary  output,  and  then  compared  to  actual 
observations.  Care  was  exercised  not  to  influence  predicted  LUF  deviations  with  prior 
knowledge  of  the  actual  event. 


Table  2  summarizes  the  X-ray  disturbances  selected  for  this  test.  A  total  of 
26  events  between  4  November  1968  and  7  August  1972  were  used.  In  the  classifi¬ 
cation  of  X-ray  bursts,  it  has  been  found  that  onlv  bursts  whose  1-8-A  flux  level  rises 
above  1  x  10  erg/cnr  •  s  produce  appreciable  SWFs  on  HF  circuits.  This  test  is 
concerned  with  two  classes  of  bursts:  small  bursts  or  Class  M  (  1  *  10  “  <  flux  <  I  - 
10  ).  and  moderate  to  large  bursts  or  class  X  ( 1  *  10-1  <  flux).  A  burst  that  peaked 

at  5  x  10~2  erg/cm^  •  s  would  be  denoted  as  an  M5  event.  Of  the  26  events.  14  were 
Class  M  and  12  were  Class  X.  The  durations  of  the  X-ray  bursts  ranged  between  15  and 
330  minutes.  There  is  an  excellent  cross  section  of  characteristics  included,  rarging  from 
very  sharp  impulsive  changes  to  very  slowly  changing  bursts.  The  smallest  burst  tested 
w'as  an  M2  and  the  largest  was  an  X4.  A  complete  cross  section  of  solar  zenith  angles  is 
represented  in  the  data  sample,  ranging  from  the  midday  case  (0  )  to  the  sunrise  and 
sunset  case,  for  which  x  is  close  to  90  .  When  the  same  event  is  tested  on  several  paths, 
the  effects  of  different  solar  zenith  angles  can  be  compared. 
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10 
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12 

13 

14 

15 

16 
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18 

19 

20 
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23 

24 

25 

26 


Date 


4 

Nov 

1969 

1  1 

Nov 

1969 

1  1 

Nov 

1969 

1  1 

Nov 

1969 

1 1 

Nov 

1969 

I  1 

Nov 

1969 

1  1 

Nov 

1969 

1  1 

Nov 

1969 

27 

Nov 

1969 

22  Jul 

1970 

12 

Aug 

1970 

16 

Nov 

1970 

16 

Nov 

1970 

16 

Nov 

1970 

17 

Nov 

1970 

17 

Nov 

1970 

17 

Nov 

1970 

17 

Nov 

1970 

17 

Nov 

1970 

18 

Nov 

1970 

2 

Aug 

1972 

2 

Aug 

1972 

2 

Aug 

1972 

4 

Aug 

1972 

7 

Aug 

1972 

7 

Aug 

1972 

Table  2.  Summary  of  events  tested. 


Size  of 

Duration 

Number  of 

Solar 

Start 

1-8- A  X-ray 

of  SWF. 

HF  Paths 

Sunspot 

Time 

Burst 

min 

Tested 

Number 

0400 

X2 

40 

1 

105 

0030 

M2 

40 

1 

105 

0130 

M7 

90 

1 

105 

0400 

XI 

40 

1 

105 

0450 

M5 

30 

1 

105 

0520 

XI 

90 

1 

105 

0710 

M5 

15 

1 

105 

0810 

X3 

30 

1 

105 

0030 

M3 

330 

1 

105 

0010 

X2 

60 

3 

103 

2000 

X3 

120 

3 

101 

0030 

X2 

90 

2 

89 

2120 

M4 

40 

2 

89 

2220 

M5 

60 

2 

89 

01  10 

M5 

50 

2 

89 

1740 

XI 

40 

2 

89 

2010 

M5 

30 

s 

89 

2120 

M2 

20 

s 

89 

2220 

M5 

40 

2 

89 

0140 

X2 

50 

2 

89 

0310 

X2 

120 

-i 

65 

1810 

M5 

70 

s 

65 

2000 

X3 

240 

s 

65 

0620 

X4 

120 

4 

65 

0200 

M3 

90 

4 

65 

0340 

M5 

70 

4 

65 

Table  3  lists  the  HF  sounder  paths  where  data  were  recorded  during  the  events  listed  in 
table  2.  Data  from  light  transmission  paths  were  tested,  spanning  the  Pacific  area  from 
California  to  Northwest  Cape,  Australia.  Here,  again,  a  variety  of  transmission  path 
characteristics  are  represented, 


Table  3.  Summary  ol  HE  transmission  paths  tested 


Path 

Path 

Length 

km 

Path 

Latitude 

Events  Tested 
(  From  Table  2  ) 

COMMSTA  Honolulu.  Hawaii, 
to  Southern  California 

4100 

W  L 
midlat 

10. 1  1.12.13.14. 

15. 16. 17.  IS.  14. 
20.2 1 .22.23 

Davis.  California,  to 

COMMSTA  Kodiak.  Alaska 

3100 

S-N 

midlat 

10. 1  1 

COMMSTA  Honolulu.  Hawaii, 
to  COMMSTA  Kodiak.  Alaska 

4100 

S-N 

midlat 

1 2. 13. 14. 15. 16. 
17.18.19.20 

COMMSTA  Guam  to 

COMMSTA  Honolulu.  Hawaii 

6000 

WE 

low  midlat 

10. 1  1.21.22.23. 
24.25.26 

COMMSTA  Guam  to 

COMMSTA  HFH.  Australia 

5300 

N-S 

transequatorial 

2. 3. 4.5.6. 7. 8. 
9.24.25.26 

COMMSTA  Guam  to 

COMMSTA  Japan 

2700 

S-N 

low  midlat 

24.25.26 

COMMSTA  Philippines  to 
COMMSTA  Japan 

3000 

S-N 

midlat 

No  tests 
performed 

COMMSTA  Philippines  to 
COMMSTA  HEH.  Australia 

4300 

N-S 

transequatorial 

1 

COMMSTA  Honolulu.  Hawaii, 
to  COMMSTA  Japan 

6000 

F.-W 

midlat 

24.25.26 

3.3  DISCUSSION  OF  ACCURACY  STUDIES 
3.3.1  Long-Range  Algorithm  Tests 

Figure  4  is  a  series  of  graphs  (a-t)  showing  predicted  LUF  (solid  lines)  and 
observed  LOF  (dashed  lines)  versus  time  for  the  26  SWF  events  resulting  from  the  solar 
X-ray  bursts  listed  in  table  3.  Each  point  plotted  in  these  graphs  represents  one  test  of 
the  SIDGRID  model. 

Figures  4a.  b.  and  c  show  three  different  SWF  situations:  a  short,  impulsive 
event  caused  by  an  X2  X-ray  burst  with  very  steep  rise  and  decay  (a);  a  very  slow,  long- 
lasting  M3  event  (b):  and  a  series  of  seven  Class  M  and  X  events  occurring  in  a  9-hour 
period  (c).  Of  interest  is  the  fact  that  these  events  occurred  in  1969.  near  the  peak  of  the 
solar  cycle,  and  the  transmission  paths  tested  are  markedly  different  from  the  path  that 
was  used  to  derive  the  original  LUF  /X-ray  flux  algorithm.  A  total  of  38  points  were 
tested  on  this  path,  with  the  error  exceeding  ±  1  MHz  in  only  one  case. 

Figures  4d  and  c  depict  two  Class  X  events  in  1970.  In  each  event,  data  from 
three  different  Eastern  Pacific  (  F.ASTPAC )  HF  transmission  paths  were  plotted.  Errors 
between  predicted  and  observed  LOF  variations  occur  in  the  decay  portion  of  the  event. 

In  figure  4d.  the  greatest  error  observed  was  on  the  Hawaii  to  California  path  between 
0100  and  0200  UT,  where  the  predicted  LUF  is  1.5  MH/  below  observed  data.  The 
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same  trend  appears  in  figure  4  for  the  same  path.  However,  in  figure  4e  for  the 
Califomia-to- Kodiak.  Alaska,  path  the  trend  is  reversed  and  the  predicted  LUF  is 
1. 5-2.0  MHz  higher  than  the  observed  LOF. 

Figures  4f  and  g  contrast  the  effects  of  a  large  X2  event  and  a  pair  of  small  Class 
M  events,  which  occurred  16  November  1970.  Figure  4f  shows  the  Hawaii-to-California 
path(W-E)  and  the  Hawaii-to- Kodiak.  Alaska,  path  (S-N).  The  decay  portion  of  the 
predicted  X3  SWF  follows  exactly  the  observed  data  on  both  paths.  On  the  M4  and  M5 
SWFs,  some  error  is  evident  and  results  from  a  lack  of  frequency  resolution  in  sounder 
data  and  also  from  a  variation  in  frequency  characteristics  of  the  predicted  LUF  algorithm 
at  the  lower  flux  level.  In  all  cases,  the  error  is  less  than  2  MHz  and.  more  tvpicallv.  it 
is  0.5  MHz. 

Figures  4h  and  i  depict  another  multiple  SWF  event  period,  which  occurred  17 
and  18  November  1970.  Again,  two  paths  are  shown:  Hawaii-to-California  and 
Hawaii-to- Alaska  events.  These  Figures  illustrate  the  different  impact  of  the  same  SWF 
on  different  transmission  paths,  depending  on  path  illumination.  At  0130  UT.  both  paths 
are  still  sunlit,  and  the  M5  X-ray  burst  produces  about  the  same  effect.  At  1800  UT.  the 
Hawaii-to-California  path  has  just  finished  sunrise,  while  the  Hawaii-to-Kodiak.  Alaska, 
path  is  still  dark.  The  XI  event  produces  a  significant  SWF  on  the  sunlit  path  and  no 
effect  on  the  darker  path.  The  SIDGRID  model  accurately  predicted  the  impact. 
Throughout  the  remaining  daylight  hours,  three  additional  Class  M  events  produced  minor 
SWFs  and  comparisons  show  good  correlation  between  predicted  and  observed  results. 
Between  0000  UT  and  0200  UT  on  18  November,  the  observed  LOF  drops  below  the 
predicted  LUF  on  both  paths.  It  should  be  noted  that  the  observed  LOF  is  anomalously 
low  for  this  time  of  day  when  compared  to  normal  daytime  LOF  trends.  In  addition, 
through  this  period,  the  1-8-A  X-ray  flux  never  dropped  below'  5  x  10“-*  erg/cm^  •  s. 
which  should  have  tended  to  keep  the  LOF  at  or  slightly  above  undisturbed  levels.  The 
final  event  occurred  during  sunset  transition  at  0230  UT.  and  again  the  model  accurately 
predicted  the  SWF  variations  over  the  two  paths.  Except  for  the  one  inconsistent  period 
between  0000  and  0200  UT.  predicted  LUFs  correlated  to  observed  LOFs  during  a 
sunrise  transition,  a  sunset  transition,  and  peak  daylight  hours.  This  demonstrates  the 
model's  capability  to  assess  SWF  impact  over  a  large  geographical  area. 

The  test  would  not  be  complete  without  a  review  of  some  of  the  SWF  events  that 
occurred  in  August  1972.  The  period  between  2  and  10  August  1972  was  the  most 
sustained  period  of  solar/ionospheric  disturbances  during  this  solar  cycle.  Figures  4 j  and 
k  show  three  SWF  disturbances  on  2  August  on  the  east-west  paths  between  Hawaii  and 
California,  and  between  Guam  and  Hawaii.  The  first  X2  event  at  0400  UT  shows  again 
the  different  SID  impact  as  a  function  of  path  longitude.  Between  Hawaii  and  California, 
the  sun  has  almost  set  and  the  impact  of  this  flare  is  minor.  However,  the  Guam-to- 
Hawaii  path  is  still  completely  illuminated,  and  the  frequency  fadeout  over  this  path 
approaches  31  MHz.  At  1900  UT.  over  this  same  path,  the  M5  X-ray  burst  tends  to 
accelerate  the  sunrise  increase  in  LOF.  w'hile  on  the  Hawaii-to-California  path,  the  M5 
event  products  a  minor  SWF.  At  2030  UT.  the  X3  event  produces  total  fadeout  over 
both  paths,  which  are  now  totally  illuminated.  The  display  used  in  the  analysis  of  these 
two  events  is  included  in  appendix  B  for  reference. 
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Figures  4  1  — o  show  the  effect  of  the  X4  X-ray  event  of  4  August  on  four  Pacific 
paths.  Note  the  difference  in  impact  over  these  paths.  The  outage  time  at  18  MHz 
varies  between  30  and  90  minutes,  and  the  model  was  able  to  accurately  predict  these 
variations  over  the  four  paths. 

Figures  4p— t  show  a  pair  of  small  Class  M  events  on  7  August,  again  for  four 
Pacific  paths.  While  the  Guam-to- Hawaii  and  Guam-to- Japan  paths  correlate 
accurately,  the  Guam-to-Japan  and  Hawaii-to-Japan  predictions  show  a  consistently 
higher  predicted  LUF  than  observed.  Of  all  the  SWF  events  and  transmission  paths 
tested,  the  plots  in  figures  4s  and  t  show  the  poorest  correlation.  There  is  one  possible 
reason  for  this  error.  First,  HF  sounder  data  received  fiom  Japan  show  a  high  daytime 
noise  environment,  making  changes  in  LOF  difficult  to  interpret.  Once  the  observed 
LOF  rises  above  16  MHz,  the  noise  contamination  no  longer  is  a  problem.  The  large 
X4  event  of  figure  4  m  shows  an  accurate  comparison  between  predicted  and  observed 
values,  while  the  smaller  M3  and  M5  events  of  figures  4s  and  t  involve  LOF  variations 
between  10  and  16  MHz,  the  part  of  the  spectrum  most  degraded  by  the  local  noise 
environment. 

Overall.  198  data-point  tests  wore  made,  comparing  SIDGR1D  predictions  to 
observed  sounder  data.  From  these  tests  it  was  determined  that: 

1.  The  mean  prediction  error  is  0.24  MHz  high. 

2.  The  error  between  and  observed  data  was  negligible  in  29%  of  the  tests. 

3.  The  probability  that  prediction  error  g  ±0.5  MHz  =48%. 

4.  The  probability  that  prediction  error  ^  ±1.0  MHz  =  83%. 

5.  The  probability  that  prediction  errors  ±2.0  MHz  =  95%. 

The  cross-sectional  data  tested  indicate  that  the  model  can  produce  valid  predictions  for 
the  entire  sunlit  hemisphere,  given  date.  time,  and  X-ray  flux.  An  accuracy  specification 
of  ±2.0  MHz  is  valid  for  X-ray  flux  levels  down  to  10“^  erg/cm^  •  s.  an  essentially 
undisturbed  condition.  For  X-ray  fluxes  of  5  x  10“^  or  greater,  the  ±0.5-MHz  error 
specification  is  more  appropriate.  When  the  initial  LUF/X-ray  flux  algorithm  was 
developed,  it  was  stipulated  that  it  was  for  a  disturbance  case  and  probably  not  valid  for 
quiet  solar  conditions  when  the  X-ray  flux  is  between  10  “4  and  10“ .  However,  these 
tests  show  that  the  SIDGRID  model  produces  usable  results  even  at  these  low  flux  levels. 

3.3.2  Short-Range  Algorithm  Tests 

The  first  performed  test  of  the  short-range  disturbance  algorithm  used  the  SWF 
disturbances  listed  in  table  1.  The  same  testing  procedures  outlined  above  were  again 
followed.  Table  4  summarizes  the  results  of  the  initial  tests.  From  the  nine  events 
tested,  the  calculated  SWF/LUF  varied  from  the  observed  SWF  over  a  range  of  -0.54 
MHz  to  ±0.20  MHz.  The  largest  deviation  occurred  in  the  12  November  1971  event, 
where  it  should  be  pointed  out  that  the  peak  X-ray  flux  was  approximated,  since  sensor 
saturation  occurred.  Also,  in  the  frequency  between  2  and  4  MHz.  the  resolution  of  the 
UPR-2  sounding  receiver  Is  0. 1  MHz.  Between  4  and  8  MHz.  the  resolution  is  0.2 
MHz.  If  these  limitations  on  the  data  are  taken  into  account,  it  is  felt  that  a  nominal 
accuracy  specification  for  the  short-range  algorithm  is  ±0.5  MHz  for  HF  transmission 
paths  of  less  than  1500  km. 
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Table  4.  Short-range  HF  LOF/ X-ray  flux  algorithm  tests.  803  km  path. 


Date 

Time 

Fx 

Observed 
SWF  LOF 

Calculated 
SWF  LUF 

Error.  1 
MH/ 

Pre-event 

LOF 

2  Jul 1971 

1845 

20.00 

0.022 

11.50 

1  1.51 

c0. 01 

5.95 

7  Ju! 1971 

1507 

62.906 

0.0063 

4.95 

4.89 

0.06 

3.675 

7  Jul 1971 

1603 

53.598 

0.0089 

6.35  | 

6.53 

■+-().  1 8 

4.75 

7  Ju! 1971 

1925 

15.0 

0.0028 

6.95 

7.02 

-t-0.07 

5.75 

7  Jul 1971 

1947 

12.192 

0.0053 

8.30 

8.31 

+0.01 

5.75 

10  Jul  1971 

1546 

54.675 

0.013 

7.15 

7.01 

-0.14  | 

3.875 

12  Sep  197! 

1624 

61.125 

0.044 

8.70 

8.31 

-0.39 

4.15 

12  Nov  1971 

1833 

65.93 

0.02 

6.55 

6.01 

-0.54 

4.75 

20  Nov  197! 

2049 

56.454 

0.018  i 

1 

7.15  | 

7.35 

+0.20  | 

357  s 

3.3.3  Complete  Algorithm  Test 

A  final  test  of  the  SIDGRID  model  was  performed  in  an  SWF  event,  for  which 
oblique  sounder  data  on  both  the  ( long-range)  Hawaii-to-La  Posta  and  ( short-range  I 
Davis-to-La  Posta  paths  were  available.  The  event  chosen  was  an  M4  X-ray  event  that 
peaked  at  00  UT  on  10  November  1971.  The  Davis-to-La  Posta  path  was  very  near 
darkness,  and  the  path  LOF"  had  already  dropped  to  2  MHz  just  prior  to  the  event.  It 
would  be  interesting  to  see  whether  the  new  short-range  algorithm  would  show  this,  along 
with  the  rise  of  the  LOF,  at  the  SWF  peak.  Figure  5  shows  excellent  agreement  between 
SIDGRID  model  predictions  and  observed  oblique  sounder  information  for  this  event. 

9-10  NOV  1971 
M4  X  RAY  SWF 


Figure  5.  Two-path  comparison  of  revised  SIDGRID 
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4.0  SUMMARY  AND  CONCLUSIONS 


1.  The  NOSC  solar  X-ray/HF  LOF  algorithms  provide  an  accurate  assessment  of 
the  change  in  HF  LOF  for  solar  X-ray  bursts  greater  than  1  -  10'  ^  erg/ cm-  •  s  and  are 
valid  over  the  entire  sunlit  hemisphere.  The  SIDGRID  model  has  the  ability  to  display 
the  disturbed  LUF  for  any  path  on  the  geographical  sunlit  hemisphere  at  any  date  and 
time.  If  the  model  is  run  with  a  continuous  date/time  clock,  the  only  required  input 
variable  is  1-8-A  X-ray  flux.  Figure  6  shows  the  SIDGRID  model  display  for  an  X-ray 
flare  starting  at  1600  UT  and  peaking  at  1615  UT.  The  display  shows  the  circuit  blacked 
out  from  about  1610  UT  to  1700  UT. 


1  MAR  1987  Hill 
10.7  CM  FLUX  ■  73 
X-RAY  FLUX  *  0. 15000 
MAGNETIC  INDEX  *2.0 
XMTR:  HA  40.00N 
ANT:  25  <0MNI  DIR) 
BEARING  < TO  SD  )  = 
RANGE  =  3825  KM  P 


mu 

■  73  SUNSPOT  NUMBER 
0. 15000 


75.00H 


=  271.4  DEG 
POUER:  1000.  W 


RCUR!  SD  32.73N  117. 17H 
ANT:  25  (OMNI  DIR) 

BEARING  <T0  HA  >  -  65.6  DEG 


FLARE  START:  16002  FLARE  PEAK:  16152 


UNIVERSAL  time 


Figure  6.  SIDGRID  model  display. 


2.  The  SIDGRID  models  HF  LUF  over  a  transmission  path  based  on  the 
following  range-dependent  parameters: 

Short-range  ^  1500  km 

Long-range  >  1500  km. 

Accuracy  of  the  long-range  algorithm  was  determined  by  comparing  SIDGRID 
predictions  to  198  observed  sounder  data  points.  Results  from  these  tests  indicate  that 
the  mean  prediction  error  was  0.24  MHz  and  the  error  between  predicted  and  observed 
data  was  negligible  in  29%  of  the  tests.  The  probability  that  the  prediction  error  was  less 
than  or  equal  to  ±1.0  MHz  was  83%  and  the  probability  of  error  less  than  +  2.0  MHz 
was  95%. 
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Accuracy  of  the  short-range  algorithm  is  shown  in  table  4,  Nine  events  were 
tested:  the  calculated  SWF  LUF  differed  from  the  observed  SWF  over  a  range  ot  u  54 
MHz  to  +0.20  MHz.  The  nominal  accuracy  specification  for  the  short-range  algorithm 
is  ±0.5  MHz  for  HF  transmission  paths  of  less  than  1500  km.  For  both  algorithms,  the 
accuracy  of  the  computed  HF  LUF  improves  as  the  solar  event  size  and  duration 
increase. 

The  sudden  ionospheric  disturbance  FORTRAN  computer  subroutine  called 
DLOF  is  listed  in  Appendix  C.  and  the  BASIC  subroutine  is  listed  in  Appendix  D. 

3.  The  SIDGRID  model  does  not  include  a  means  to  adjust  the  resultant  l.L'Fs 
for  system  parameters  such  as  transmitted  power,  antenna  gains,  required  signal -to  noise 
ratio,  and  background  radio  noise.  Under  the  full  impact  of  an  SWF.  the  main  contribu 
tion  to  the  LUF  is  the  D-region  absorption.  However,  as  the  D- region  begins  to  rccovei 
from  the  effects  of  the  SWF.  the  transition  to  quiet  time  LUF  prediction  would  be 
smoother  if  methods  were  developed  to  include  these  effects. 

4.  The  SIDGRID  model  is  strictly  a  LUF  prediction  model.  It  does  not  calculate 
the  D-region  absorption  resulting  from  an  SWF.  Such  a  model  is  desirable  so  that  the 
decrease  in  available  signal  strength  between  the  LUF  and  MUF  resulting  from  an  SWT 
can  be  determined. 

5.  The  SIDGRID  model  is  a  two-range  model  (short  range  <  1500  km.  long 
range  >  1500  km).  An  attempt  should  be  made  to  keep  the  model  continuous  as  a 
function  of  range. 
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APPENDIX  A: 
THEORETICAL  BASIS 


It  is  appropriate  to  look  briefly  at  why  the  simple  logarithm  used  to  produce  the 
SIDGRID  works  as  well  as  it  does.  In  references  10  and  I  I .  it  was  determined  that  the 
undisturbed  LUF  is  approximated  by: 

LUF  =  (Constant)  (cos  \V4  . 

where  Constant  =  a  function  of  D-region  nondeviative  absorption  (  A)  plus  transmitted 
power,  system  signal-to-noisc  ratio,  launch  angle,  and  antenna  patterns. 

In  reference  10.  it  also  was  stated  that  the  absorption  (  A)  could  be  expressed  by: 


where 

N(C|  ~  electron  density  of  the  D  region 
t  frequency  of  operation. 

The  D  region  electron  density  can  be  defined  as  consisting  of: 

Nfc)  =  N<  e )  X  +  N(e)UV 


where 

N.  ,  v  electron  density  due  to  X-rax  emission 

U  I  A 

N(e)yv  electron  density  due  to  ultraviolet  Lyman  Alpha  (  12  lb  -  A)  emission. 

N(e)UV  remains  essentially  constant  during  quiet  and  disturbed  solar  conditions. 

N(  e )  x  's  a  secondary  contributor  to  absorption  during  undisturbed  periods  and  becomes 
the  primary  contributor  during  disturbed  periods  The  SIDGRID  algorithm  was 
empirically  developed  on  the  basis  of  observed  changes  in  X-ray  emission  and  HE  l  Of 
Because  this  relationship  was  established  for  the  cases  in  which  solar  X  ray  emissions  arc 
the  principal  cause  of  HE  absorption,  the  accuracy  of  the  algorithm  should  increase 
proportionally  to  the  si/e  of  the  X  ray  burst  In  other  words,  the  SIDGRID  model  is 
predicated  on  the  expression: 

LUF  •  I  Absorp  \y  r;n  )  '  (cosy)4 

As  the  magnitude  of  the  input  X  ray  emissions  decreases,  the  HI  system  parameters 
increase  in  importance,  complicating  the  constant  calculation,  and  tor  the  undisturbed 
case  the  HF  LOE  is  expressed  by 

LUF  ■  (HE  system  path  parameters  •  Absorp  j  y  I  ■  i  cos  yi'* 

Therefore,  for  extremely  low  levc  Is  of  X  ray  emission,  it  can  he  expected  that  the 
accuracy  of  the  SIDGRID  model  will  Jinumsh  somewh.it 


APPENDIX  B: 

TEST  OF  DOUBLE  SWF  ON  2  AUGUST  1972 


Figure  B-l  shows  a  composite  ot  the  HE  propagation  between  Hawaii  and  Guam 
for  the  24-hour  period  of  2-3  August  1972.  the  1  8  A  X  ra\  hursts  that  caused  the  two 
SWFs  and  (shaded  area)  the  predicted  SWF  derived  from  the  SIDGRID  model  lime 
resolution  is:  one  line  equals  a  10-minute  average  of  the  data.  Signal  amplitude  is  25  d 
above  the  receiver  threshold. 


Figure  B-1.  24-hour  propagation  characteristics  Guam  to  Hawaii  2 -3  August  19  7?  Time 
resolution-  1  line  -  10  mm  average  Signal  amplitude-  \  25  dB  above  receiver  threshold 


APPENDIX  C: 

DLOF  —  FORTRAN  SUDDEN  IONOSPHERIC 
DISTURBANCE  SUBROUTINE 


subroutine  dlof  (  tlat,  tlon,  rlat,  rlon,  .xflux,  cpnt,  lof  ) 
cp 

c  subroutine  dlof 
c 

c  call  dlof(tlat, tlon, rlat, rlon, xflux, cpnt, lof) 

c 

c  this  routine  computes  the  luf  for  solar  disturbed  conditions  (i.e. 

c  xray  flux  >  1 .8 3e - 3 .  tlat  'tlon  &  rlat/rlon  are  the  end  point 

c  coordinates  in  radians,  xflux  is  the  current  xray  flux,  cpnt  is  an 

c  eight  element  array  containing  the  path  length  and  coordinates  of 

c  control  points  along  the  path  (see  subroutine  path  for  complete 

c  description),  this  routine  returns  lof  which  is  the  luf  in  mhz. 

c  this  routine  assumes  west  longitudes, 
c 

c  subroutines  and  functions  used:  minpt 

c  newton 

c 

c  common  blocks:  sun 

cz 

real  cpnt(8),  lof 

c 

common  sun  slat,  slon 

c 

lof  =  0.5 

call  minpt(  rlat,  rlon,  cpnt.  clat,  cion,  chi  ) 

dist  =  cpnt(  1  )*637 1 .0 

if  (  chi  gt.  1 .57  )  go  to  100 

flux  =  xflux 

c 

if  (  dist  It.  3500.0  )  go  to  1 1 
call  newtont  chi,  flux,  lof  ) 
go  to  100 

11  x  =  sqrt(  (  flux*cos(  chi  )**3  ),'!.03856e-6  ) 

20  lof  =  sqrt(  x  ) 

theta  =  (  dist/6371.0  )*0.5 

alfa  =  (  cos(  theta  )  -  0.96224  )/sin(  theta  ) 

alfa  =  atan(  alfa  ) 

alfa  =  acos  (  0.989  l*cos(  alfa  )  ) 

lof  =  lof*sqrt(  0.5368/ sin(  alfa  )  ) 

100  continue 

lof  =  aminl(  amaxl(  lof,  0.5  )  ,  50.0  ) 

return 

end 


APPENDIX  D: 

SIDGRID—  BASIC  SUDDEN  IONOSPHERIC 
DISTURBANCE  SUBROUTINE 


4750  SUB  Sidgrid 

4760  A=S0 

4770  B=S1 

4780  C=Z3 

4790  D=Z4 

4800  CALL  Gcraz 

4810  C7=R 

4820  C=Z5 

4830  D=Z6 

4840  CALL  Gcraz 

4850  C7=C7  MIN  R 

4860  IF  C7<P0  THEN  4890 

4870  L=2 

4880  GO  TO  5320 

4890  FOR  1=1  TO  9 

4900  C=Z0(2*I+2) 

4910  IF  C=-999  THEN  4960 

4920  D=Z0(2*I+3) 

4930  CALL  Gcraz 

4940  C7=C7  MIN  R 

4950  NEXT  I 
4960  C7=COS(C7) 

4970  D=Z0(I)*R2 

4980  IF  D=>3500  THEN  5070 

4990  X=SQR(SQR(  F4*C7' 3  1 .038561.-6 1 ) 

5000  T=D  12742 

5010  A  =  ATN((COSn  )-0.96224i  SIN*  ID 
5020  A  =0.9891  *COS(  A ) 

5030  A  =  A  MAX  -1  MIN  I 
5040  A=ACS(  A ) 

5050  L=X*SQR(0.5368  S1N(A)) 

5060  GO  TO  5310 
5070  B=  I  C7" 2 

5080  X0=(F4  0.1038+150)  *  10+1!) 
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C=  1  .OE-4 

X=XO 

E=X 

CALL  Newton 

F=100*C 

FOR  1=1  TO  20 

IF  G=0  THEN  5300 
IF  H=0  THEN  5290 
X 1 =G/H 
X=X-X1 
E=X 

CALL  Newton 
E=C 

R=ABS(X) 

IF  R<=  1  THEN  5250 
E=E*R 

IF  ABS(X  1  )>E  THEN  5270 
IF  ABS(G)<=F  THEN  5300 
NEXT  I 
GO  TO  5300 
X=32 
L=X 

CALL  Adjluf 
END  SUB 
SUB  Newton 
P=(1+B/10)*E 
Q=0.849*(P-  1  5.6) 

Q=Q-INT(Q/P1  )*PI 
G=0.01038*(P- 1 5)-0.003*SIN(Q) 
G=G-F4 

H=(0.0 1038-0. 0025473*COS(Q))*P  E 
END  SUB 


